Quark-Gluon Bags with Surface Tension 
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The temperature and chemical potential dependent surface tension of bags is introduced into the 
gas of quark-gluon bags model. This resolves a long standing problem of a unified description of the 
first and second order phase transition with the cross-over. Such an approach is necessary to model 
the complicated properties of quark-gluon plasma and hadronic matter from the first principles 
of statistical mechanics. The suggested model has an exact analytical solution and allows one to 
rigorously study the vicinity of the critical endpoint of the deconfinement phase transition. The 
existence of higher order phase transitions at the critical endpoint is discussed. In addition, we 
found that at the curve of a zero surface tension coefficient there must exist the surface induced 
phase tranition of the 2 nd or higher order, which separates the pure quark gluon plasma (QGP) 
from the cross-over states, that are the mixed states of hadrons and QGP bags. Thus, the present 
model predicts that the critical endpoint of quantum chromodynamics is the tricritical endpoint. 
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I. INTRODUCTION 

Investigation of the strongly interacting matter proper- 
ties observed in relativistic nuclear collisions has reached 
the stage when the predictions of the lattice quantum 
chromodynamics (QCD) can be checked experimentally 
on the existing data and future mesurements at BNL 
RHIC, CERN SPS, and GSI FAIR. However, a compari- 
son of the theoretical results with the experimental data 
is not straightforward because during the collision pro- 
cess the matter can have several phase transformations 
which are difficult to model. The latter reason stimu- 
lated the development of a wide range of phenomenologi- 
cal models of the strongly interacting matter equation of 
state which are used in dynamical simulations. 

One of these models is the famous bag model [1] which 
treats the hadrons as the bags of QGP confined inside a 
hadron with help of bag pressure. The bag model is able 
to simultaneously describe the hadron mass spectrum, 
i.e. the hadron masses and their proper volumes, and 
the properties of the deconfined phase 0). This success 
led to a development of the statistical model of QGP, 
the gas of bags model (GBM) which itself con- 

tains two well-known models of deconfined and confined 
phases: the bag model of QGP 2] and the hadron gas 
model Q. There were hopes 0] that an exact analyt- 
ical solution of the GBM found in Q could be helpful 
in understanding the properties of strongly interacting 
matter. However, this solution does not allow one to in- 
troduce the critical end point of the strongly interacting 
matter phase diagram. Also, a complicated construc- 
tion of the line, along which the phase transition order 
gradually increases, suggested in [3], does look too arti- 
ficial. Therefore, the present GBM formulation lacks an 
important physical input and is interesting only as a toy 
example which can be solved analytically. 

On the other hand, the models, which can correctly 
reproduce the expectation [g, Qj, [lCj that the end point of 
the I s * order phase transition (PT) line to QGP should 
be the 2 nd order critical point, are indeed necessary for 



heavy ion phenomenology. In addition, such phenomeno- 
logical models can provide us with the information about 
the phase structure and equation of state of strongly in- 
teracting matter which is located between the critical 
endpoint and the region of the color superconductivity 
because such an information is unavailable otherwise. 
Therefore, the present work is devoted to the extension 
of the GMB. We think that the GMB can be drastically 
improved by the inclusion of such a vitally important 
element as the surface tension of the quark-gluon bags. 

The dynamical surface tension of the quark-gluon bags 
was estimated long ago [ll| , but it was never used in sta- 
tistical description of the equation of state. Moreover, 
the estimate of the bag surface tension made in [ll[ is 
negligible for u and d quarks of and, hence, can be safely 
neglected in our treatment. Furthermore, although the 
influence of the surface energy of the QGP bag proper- 
ties in the vacuum was discussed recently [l2[ , the surface 
tension of large bags was not included into statistical de- 
scription of QGP. Therefore, the present paper is devoted 
to the investigation and analysis of the critical properties 
of the model of quark-gluon bags with surface tension 
(QGBST model hereafter). 

In statistical mechanics there are several exactly solv- 
able cluster models with the 1 st order PT which describe 
the critical point properties very well. These models arc 
built on the assumptions that the difference of the bulk 
part (or the volume dependent part) of free energy of two 
phases disappears at phase equilibrium and that, in ad- 
dition, the difference of the surface part (or the surface 
tension) of free energy vanishes at the critical point. The 
most famous of them is the Fisher droplet model (FDM) 
[HI, 03] which has been successfully used to analyze the 
condensation of a gaseous phase (droplets of all sizes) 
into a liquid. The FDM has been applied to many differ- 
ent systems, including nuclear multifragmentation [l5|, 
nucleation of real fluids [16[ , the compressibility factor of 
real fluids fl7l , bags of the Ising model [HI and percola- 
tion bags 
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On the basis of the statistical multifragmentation 
model (SMM) [2(| commonly used to study nuclear mul- 
tifragmentation, there was recently formulated a simpli- 
fied SMM version which was solved analytically both 
for infinite [U M, and for finite [H, [H, [25| volumes 
of the system. In the SMM the surface tension tem- 
perature dependence differs from that one of the FDM, 
but it was shown [22| that the value of Fisher exponent 
tsmm = 1-825 ± 0.025, which contradicts to the FDM 
value tfdm ~ 2.16, is consistent with ISiS Collaboration 
data [26[ and EOS Collaboration data (27} • Lately, our 
analytical results [22| were confirmed by the numerical 
studies [Hil] . 

Such an experimentally obtained range of the r index 
is of a principal importance because it gives a very strong 
evidence that the SMM, and, thus, the nuclear matter, 
has a tricritical endpoint rather than a critical endpoint 

MM. 



This success of the SMM initiated the studies of the 
surface partitions of large clusters within the Hills and 
Dales Model [3(| HH and led to a discovery of the origin 
of the temperature independent surface entropy similar 
to the FDM. As a consequence, the surface tension co- 
efficient of large clusters consisting of the discrete con- 
stituents should linearly depend on the temperature of 
the system 30] and must vanish at the critical endpoint. 
However, the present formulation of the Hills and Dales 
Model [20, [3l|, which successfully estimates the upper 
and lower bounds of the surface deformations of the dis- 
crete physical clusters, does not look suitable for quark- 
gluon bags. Therefore, in this work we assume a certain 
dependence of the surface tension coefficient on tempera- 
ture and baryonic chemical potential, and concentrate on 
the impact of surface tension of the quark-gluon bags on 
the properties of the deconfinement phase diagram and 
the QCD critical endpoint. A discussion of the origin of 
the surface tension is a subject of our future work. 

Here we will show that the existence of a cross-over at 
low values of the baryonic chemical potential along with 
the 1 st order deconfinement PT at high baryonic chemical 
potentials leads to the existence of an additional PT of 
the 2 nd or higher order along the curve where the surface 
tension coefficient vanishes. Thus, it turns out that the 
QGBST model predicts the existence of the tricritical 
rather than critical endpoint. 

The paper is organized as follows. Sect. II contains 
the formulation of the basic ingredients of the GBM. In 
Sect. Ill we formulate the QGBST model and analyze all 
possible singularities of its isobaric partition for vanishing 
baryonic densities. This analysis is generalized to non- 
zero baryonic densities in Sect. IV. Sect. V is devoted 
to the analysis of the surface tension induced PT which 
exists above the deconfinement PT. The conclusions and 
research perspectives are summarized in Sect. V. 



II. BASIC INGREDIENTS OF THE GBM 

To remind the basic ingredients of the GBM let us 
consider the Van der Waals gas consisting of n hadronic 
species, which are called bags in what follows, at zero 
baryonic chemical potential. Its grand canonical parti- 
tion (GCP) is given by 
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is the particle density of bags of mass m& and eigen vol- 
ume Vk and degeneracy gt- Using the standard technique 
of the Laplace transformation [JHH with respect to vol- 
ume, one obtains the isobaric partition: 



Z(s,T) = / dVexp{-sV) Z(V,T) = 



(2) 



[s-F(s,T)} 
o 

n 

with F(s, T) = exp (— Vjs) gj<fi(T,mj) . (3) 
i=i 

From the definition of pressure in the grand canonical 
ensemble it follows that, in the thermodynamic limit, the 
GCP of the system behaves as Z(V, T) ~ exp [pV/T]. An 
exponentially increasing Z(V, T) generates the rightmost 
singularity s* = p/T of the function Z(s,T) in variable 
s. This is because the integral over V in Eq. diverges 
at its upper limit for s < p/T. Therefore, the rightmost 
singularity s* of Z{s,T) gives us the system pressure: 



P(T) 



T lim 

V^oo 



In Z(V, T) 



V 



= T s*{T) 



(4) 



The singularity s* of Z(s,T) 
the transcendental equation 



can be calculated from 
|2jl: 



s*(T) = F(s*,T) 



(5) 



As long as the number of bags, n, is finite, the only pos- 
sible singularities of Z(s,T) ^ are simple poles. For 
example, for the ideal gas [n = l\Vi = in Eq. ([5])) 
s* = gi<j)(T,mi) and thus from Eq. (g]) one gets p — 
Tgi4>(T,m{) which corresponds to the grand canonical 
ensemble ideal gas equation of state for the particles of 
mass mi and degeneracy g\. 

However, in the case of an infinite number of sorts of 
bags an essential singularity of Z(s, T) may appear. This 
property is used in the GBM: to the finite sum over dif- 
ferent bag states in {2]) the integral Jf? dmdv...p(m,v) 
is added with the bag mass- volume spectrum, p(m,v), 
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which defines the number of bag states in the mass- 
volume region [m,v;m + dm,v + dv]. In this case the 
function F(s, T) in Eqs. ^ and ([5]) should be replaced 
by 

71 

F(s,T) = F H (s,T) + F Q (s,T) = ^ fte -^^(T,m,) 

oo oo 



+ J dv J dm p(m,v) exp(—sv)4>(T,m) . (6) 

V M a +Bv 

The first term of Eq. ©, Fh, represents the contri- 
bution of a finite number of low-lying hadron states. 
This function has no s-singularities at any temperature 
T and can generate a simple pole of the isobaric par- 
tition, whereas the mass-volume spectrum of the bags 
Fq(s, T) can be chosen to generate an essential singular- 
ity sq(T) = pq(T)/T which defines the QGP pressure 
Pq (T) at zero baryonic densities [1, H2] • 

The mass-volume spectrum is the generalization of 
the exp onential mass spectrum introduced by Hagedorn 
[33l Hj] . The usage of the grand canonical description for 
the exponential mass spectrum was recently strongly crit- 
icized [H, |H, H3] because of the thermostatic properties 
of this spectrum. Fortunately, the Van der Waals repul- 
sion compensates the growing part of the mass-volume 
spectrum and, hence, the criticism of Refs. [HI, H^, H3] 
is irrelevant to the present model. 

There are several possibilities to parameterize the 
mass- volume spectrum p(m,v). Thus, in the simplest 
case one can assume that for heavy resonances their mass 
and eigen volume are proportional, i.e. the spectrum 
p(m, v) contains the function Sim — v Const). An alter- 
native choice was suggested in [32| , but in either case the 
resulting expression for the continuum spectrum of the 
GBM F Q (s,T) can be cast as 



F Q (s,T) = u{T) 



dv eM-v(s-s Q (T))} ^ (?) 
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where u(T), r > are the model parameters. The QGP 
pressure p(T) = Tsq(T) can be parameterized in many 
ways. For instance, the MIT bag model equation of state 
[H corresponds [H to s Q (T) = § a Q T 3 - f and u(T) = 

C*7t- 1 ( 7q +1/2 T 4 + 45 (<j q T 4 + B) 3/2 . Here B denotes the 

bag constant, ctq = |j y is the Stefan-Boltzmann 
constant counting gluons (spin, color) and (anti-)quarks 
(spin, color and u, d, s-flavor) degrees of freedom; and 
the constants C, 5 < 0, V ss 1 fm 3 and M a k, 2 GeV are 
the parameters of the mass-volume spectrum. 

III. THE ROLE OF SURFACE TENSION 

At the moment the particular choice of function 
Fq(s,T) ((7|) is not important. The key point for our 
study is that it should have the form of Eq. ([7]) which has 



a singularity at s = sq because for s < sq the integral 
over dv diverges at its upper limit. Note that the expo- 
nential in |(7J) is nothing else, but a difference of the bulk 
free energy of a bag of volume v, i.e. —Tsv, which is un- 
der external pressure Ts, and the bulk free energy of the 
same bag filled with QGP, i.e. —Tsqv. At phase equi- 
librium this difference of the bulk free energies vanishes. 
Despite all positive features, Eq. lacks the surface 
part of free energy of bags, which will be called a surface 
energy hereafter. In addition to the difference of the bulk 
free energies the realistic statistical models which demon- 
strated their validity, the FDM [H[ and SMM [20J , have 
the contribution of the surface energy which plays an 
important role in defining the phase diagram structure 
[2lL [25j . Therefore, we modify Eq. by introducing 
the surface energy of the bags in a general fashion [22j : 

00 

Fq=u{t) r dv ^m^h^m^i, (8) 
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where the ratio of the temperature dependent surface 
tension coefficient to T (the reduced surface tension co- 
efficient hereafter) which has the form a(T) — Zf- ■ 

2k+l 

(k — 0,1,2,...). Here a > can be a 

smooth function of the temperature, but for simplicity 
we fix it to be a constant. For k = the two terms in 
the surface (free) energy of a u-volume bag have a simple 
interpretation [1 31 ] : thus, the surface energy of such a bag 
is oqv* , whereas the free energy, which comes from the 



surface entropy <J T a 



is 



- cep 



Note that the 



surface entropy of a v- volume bag counts its degeneracy 
factor or the number of ways to make such a bag with all 
possible surfaces. This interpretation can be extended to 
k > on the basis of the Hills and Dales Model [13, HH- 
In choosing such a simple surface energy parameteri- 
zation we follow the original Fisher idea [13| which al- 
lows one to account for the surface energy by considering 
some mean bag of volume v and surface . The consid- 
eration of the general mass- volume-surface bag spectrum 
we leave for the future investigation. The power x < 1 
which describes the bag's effective surface is a constant 
which, in principle, can differ from the typical FDM and 



SMM value 



This is so because near the deconfinc- 



ment PT region QGP has low density and, hence, like in 



38 1 , the non-sperical bags 
38|) can be favorable. A 



the low density nuclear matter 
(spaghetti-like or lasagna-like 
similar idea of "polymerization" of gluonic quasiparticles 
was introduced recently (39j |. 

The second essential difference with the FDM and 
SMM surface tension parameterization is that we do not 
require the vanishing of a(T) above the CEP. As will 
be shown later, this is the most important assumption 
which, in contrast to the GBM, allows one to naturally 
describe the cross-over from hadron gas to QGP. Note 
that negative value of the reduced surface tension co- 
efficient cr(T) above the CEP does not mean anything 
wrong. As we discussed above, the surface tension coef- 
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ficient consists of energy and entropy parts which have 



opposite signs [13J, [30j, [31|. Therefore, a(T) < does 



not mean that the surface energy changes the sign, but 
it rather means that the surface entropy, i.e. the loga- 
rithm of the degeneracy of bags of a fixed volume, simply 
exceeds their surface energy. In other words, the num- 
ber of non-spherical bags of a fixed volume becomes so 
large that the Boltzmann exponent, which accounts for 
the energy "costs" of these bags, cannot suppress them 
anymore. 

Finally, the third essential difference with the FDM 
and SMM is that we assume that the surface tension 
in the QGBST model happens at some line in fis — T 
plane, i.e. T cep = T cep (/is)- However, in the subsequent 
sections we will consider T cep = Const for simplicity, and 
in Sect. V we will discuss the necessary modifications of 
the model with T cep = T cep (/is)- 

The surface energy should, in principle, be introduced 
into a discrete part of the mass- volume spectrum Fh , but 
a successful fitting of the particle yield ratios [|| with the 
experimentally determined hadronic spectrum Fh does 
not indicate such a necessity. 

According to the general theorem [3J] the analysis of 
PT existence of the GCP is now reduced to the analy- 
sis of the rightmost singularity of the isobaric partition 
([2]) . Depending on the sign of the reduced surface tension 
coefficient, there are three possibilities. 

(I) The first possibility corresponds to <x(T) > 0. Its 
treatment is very similar to the GBM choice ([7]) with 
t > 2 3]. In this case at low temperatures the QGP 
pressure Tsq(T) is negative and, therefore, the right- 
most singularity is a simple pole of the isobaric partition 
s* = s H {T) = F(s H (T),T) > s Q (T), which is mainly 
defined by a discrete part of the mass-volume spectrum 
Fh(s,T). The last inequality provides the convergence 
of the volume integral in ([5]) (see Fig. 1). On the other 
hand at very high T the QGP pressure dominates and, 
hence, the rightmost singularity is the essential singu- 
larity of the isobaric partition s* = sq(T). The phase 
transition occurs, when the singularities coincide: 



frp \ _ Ph(T c ) . 

Sh[T c ) = — — = s Q (T c 



pq{ t c) 



(9) 



which is nothing else, but the Gibbs criterion. The graph- 
ical solution of Eq. ([5]) for all these possibilities is shown 
in Fig. 1. Like in the GBM [3,0], the necessary condition 
for the PT existence is the finiteness of Fq(sq(T), T) at 
s = sq (T) . It can be shown that the sufficient conditions 
are the following inequalities: Fq(sq(T),T) > sq(T) for 
low temperatures and F(sq(T), T) < sq(T) for T — > oo. 
These conditions provide that at low T the rightmost sin- 
gularity of the isobaric partition is a simple pole, whereas 
for hight T the essential singularity sq(T) becomes its 
rightmost one (see Fig. 1 and a detailed analysis of case 
Ms + 0). 

The PT order can be found from the T-derivatives of 




FIG. 1: Graphical solution of Eq. ([5]) which corresponds to a 
PT. The solution of Eq. ([5]) is shown by a filled hexagon. The 
function F(s,£) is shown by a solid curve for a few values of 
the parameter £. The function F(s,£) diverges for s < sq(£) 
(shown by dashed lines), but is finite at s = sq(£) . At 
low values of the parameter £ = £a, which can be either T 
or fiB, the simple pole sh is the rightmost singularity and 
it corresponds to hadronic phase. For £ = £b 2> £a the 
rightmost singularity is an essential singularity s — Sq(£b), 
which describes QGP. At intermediate value £ = £c both 
singularities coincide sh(£c) = s q(£c) and this condition is 
a Gibbs criterion. 



sh(T). Thus, differentiating ([5]) one finds 
G + u/C r _i(A, — a) ■ s'q 



(10) 



1 + ulCr-i(A,~a) 
where the functions G and /C r _ a (A, —a) are defined as 



G = F' H + ~^Fq + (T 7^- T T\ a T T) "Er-„(A, -a) , (11) 



IT, 



/C r - Q (A, -a) 



dv 



exp [-Av - a(T)v x ] 



(12) 



where A = Sh — sq. 

Now it is easy to see that the transition is of the 1 st 
order, i.e. s'q(T c ) > s' H (T c ), provided er(T) > for any r. 
The 2 nd or higher order phase transition takes place pro- 
vided Sq(T c ) = s H (T c ) at T = T c . The latter condition is 
satisfied when /C T _i diverges to infinity at T — > [T c — 0), 
i.e. for T approaching T c from below. Like for the 
GBM choice ([7]), such a situation can exist for cr{T c ) = 
and | < t < 2. Studying the higher T-derivatives of 
sh(T) at T c , one can show that for a(T) = and for 
(n + l)/n < r < nj (n — 1) (n — 3, 4, 5, ...) there is a n th 
order phase transition 

s H {T c ) = s Q (T c ) , s' H (T c ) = s' Q (T c ) , ... 

■^{Tc) - s%- l \T c ) , s%\T c ) + S W(T C ) , (13) 

with s^iTc) = oo for (n + l)/n < t < n/(n — 1) and 
with a finite value of sff 1 (T c ) for t = (n + l)/n. 
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FIG. 2: Graphical solution of Eq. ([5} which corresponds 
to a cross-over. The notations are the same as in Fig. 1. 
Now the function F(s,£) diverges at s = sq(£) (shown by 
dashed lines) . In this case the simple pole sh is the rightmost 
singularity for any value of £. 



(II) The second possibility, <j(T) = 0, described in 
the preceding paragraph, does not give anything new 
compared to the GBM [3j, |7f. If the PT exists, then the 
graphical picture of singularities is basically similar to 
Fig. 1. The only difference is that, depending on the PT 
order, the derivatives of F(s, T) function with respect to 
s should diverge at s — sq(T c ). 

(III) A principally new possibility exists for T > T cep , 
where a{T) < 0. In this case there exists a cross-over. 
Its existence can be shown as follows. Let us solve the 
equation for singularities ([5]) graphically (see Fig. 2). For 
a(T) < the function F Q (s,T) diverges at s = sq(T). 
On the other hand, the partial derivatives ^Es^il} < o 

and < are always negative. Therefore, the 

function F(s,T) = F H (s,T) + Fq(s,T) is a monotoni- 
cally decreasing function of s, which vanishes at s — > oo. 
Since the left hand side of Eq. ([5|) is a monotonically in- 
creasing function of s, then there can exist a single inter- 
section s* of s and F(s, T) functions. Moreover, for finite 
sq (T) values this intersection can occur on the right hand 
side of the point s — sq(T), i.e. s* > sq(T) (see Fig. 
2). Thus, in this case the essential singularity s — sq(T) 
can become the rightmost one for infinite temperature 
only. In other words, the pressure of the pure QGP can 
be reached at infinite T, whereas for finite T the hadronic 
mass spectrum gives a non-zero contribution into all ther- 
modynamic functions. 

Note also that all these nice properties would vanish, 
if the reduced surface tension coefficient is zero or posi- 



tive above T cep . This is one of the crucial points of the 
present model which puts forward certain doubts about 
the vanishing of the reduced surface tension coefficient in 
the FDM pjf and SMM (20|. These doubts are also sup- 



ported by the first principle results obtained by the Hills 
and Dales Model [30, ElO], because the surface entropy 
simply counts the degeneracy of a cluster of a fixed vol- 
ume and it does not physically affect the surface energy 
of this cluster. 

IV. GENERALIZATION TO NON-ZERO 
BARYONIC DENSITIES 

The possibilities (I)- (III) discussed in the preceding 
section remain unchanged for non-zero baryonic num- 
bers. The latter should be included into consideration 
to make our model more realistic. To keep the presen- 
tation simple, we do not account for strangeness. The 
inclusion of the baryonic charge of the quark-gluon bags 
does not change the two types of singularities of the iso- 
baric partition @ and the corresponding equation for 
them ([5]), but it leads to the following modifications of 
the Fjj and Fq functions: 



3=1 

Fq(s,T,h b ) = u(T,n B ) 



x J dv 



exp [(sq(T, fi B ) — s)v — a^v" 



(14) 



(15) 



Here the baryonic chemical potential is denoted as /j,b, 
the baryonic charge of the j-th hadron in the discrete part 
of the spectrum is bj . The continuous part of the spec- 
trum, Fq can be obtained from some spectrum p(m, v, b) 
in the spirit of Ref . [321 ] , but this will lead us away from 
the main subject. 

The QGP pressure pq = Tsq(T, (1b) can be also cho- 
sen in several ways. Here we use the bag model pressure 



95 10 /mb\ 2 

9 + n 2 { rp j 



{ T ) + 9tt 4 ( T ) 



B, (16) 



but the more complicated model pressures, even with the 
PT of other kind like the transition between the color 
superconducting QGP and the usual QGP, can be, in 
principle, used. 

The sufficient conditions for a PT existence are 

F(s Q (T,ti B =Q),T,(i B = Q)>s Q (T,fi B =Q), (17) 
F(sq(T,/j, b ),T,hb) < s q (T,hb) , for all fj, > fi A - (18) 



The condition (I17| provides that the simple pole singular- 
ity s* = sh(T, (is — 0) is the rightmost one at vanishing 
fj, B = and given T, whereas the condition (|18|) ensures 
that s* = SQ(T,fis) is the rightmost singularity of the 
isobaric partition for all values of the baryonic chemical 
potential above some positive constant \xa- This can be 
seen in Fig. 1 for [Ib being a variable. Since F(s, T,fi B ), 
where it exists, is a continuous function of its parame- 
ters, one concludes that, if the conditions (fTT]) and (fT8]). 
are fulfilled, then at some chemical potential A*g(T) the 
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both singularities should be equal. Thus, one arrives at 
the Gibbs criterion but for two variables 



sh (T, (i c B (T) ) = sq (T, p c B (T) ) 



(19) 



It is easy to see that the inequalities (fTT)) and (fT5|) 
are the sufficient conditions of a PT existence for more 
complicated functional dependencies of Fh{s, T, fis) and 
Fq(s,T, fis) than the ones used here. 
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FIG. 3: A schematic picture of the deconfinement phase 
transition diagram in the plane of baryonic density pB and T 
for the 2 nd order PT at the critical endpoint (CEP), i.e. for 
| < t < 2. For the 3 rd (or higher) order PT the boundary 
of the mixed and hadronic phases (dashed curve) should have 
the same slope as the boundary of the mixed phase and QGP 
(solid curve) at the CEP. 



ture of the phase diagram in the variables temperature- 
density has, to our knowledge, never been observed, we 
conclude that the case r > 2 is unrealistic (see Fig. 4). 
Note that a similar phase diagram exists in the FDM 
with the only difference that the boundary of the mixed 
and liquid phases (the latter in the QGBST model cor- 
responds to QGP) is moved to infinite particle density. 



CELine 




Pb 

FIG. 4: A schematic picture of the deconfinement phase 
transition diagram in the plane of baryonic density ps and 
T for t > 2. The critical endpoint in the ps — T plane 
generates the critical end line (CELine) in the pB — T plane 
shown by the thick horizontal line. This occurs because of the 
discontinuity of the partial derivatives of sh and sq functions 
with respect to ps and T. 



For our choice ([11]), flTSJ) and {TBJ of F H (s,T, fi B ) and 
Fq(s,T, fis) functions the PT exists at T < T cep , be- 
cause the sufficient conditions ([T71) and (fT8|) can be eas- 
ily fulfilled by a proper choice of the bag constant B and 
the function m(T, fig) > for the interval T < T up with 
the constant T up > T cep . Clearly, this is the I s * order 
PT, since the surface tension is finite and it provides the 
convergence of the integrals ifTTj) and (fl"2")) in the expres- 
sion pTJj) . where the usual T-derivatives should be now 
understood as the partial ones for ps = const. 

Assuming that the conditions (|17[) and (fl"5|) arc fulfilled 
by the correct choice of the model parameters B and 
u(T,fiB) > 0, one can see now that at T = T cep there 
exists a PT as well, but its order is defined by the value 
of t. As was discussed in the preceding section for | < 
t < 2 there exists the 2 nd order PT. For 1 < r < § there 
exist the PT of higher order, defined by the conditions 
formulated in Eq. (|13|) . This is a new possibility, which, 
to our best knowledge, does not contradict to any general 
physical principle (see Fig. 3). 

The case r > 2 can be ruled out because there must ex- 
ist the first order PT for T > T cep , whereas for T < T cep 
there exists the cross-over. Thus, the critical endpoint in 
T — fiB plane will correspond to the critical interval in the 
temperature-baryonic density plane. Since such a struc- 



V. SURFACE TENSION INDUCED PHASE 
TRANSITION 

Using our results for the case (III) of the preceding 
section, we conclude that above T cep there is a cross- 
over, i.e. the QGP and hadrons coexist together up to 
the infinite values of T and/or ps- Now, however, it 
is necessary to answer the question: How can the two 
different sets of singularities that exist on two sides of 
the line T = T cep provide the continuity of the solution 
ofEq. ©? 

It is easy to answer this question for fis < Ms(^cep) be- 
cause in this case all partial T derivatives of sh{T, pb), 
which is the rightmost singularity, exist and are finite 
at any point of the line T = T cep . This can be seen 
from the fact that for the considered region of parameters 
sj^(T, fis) is the rightmost singularity and, consequently, 
sh(T,pb) > sq(T,plb)- The latter inequality provides 
the existence and finiteness of the volume integral in 
Fq(s,T, ps)- In combination with the power T depen- 
dence of the reduced surface tension coefficient <r(T) the 
same inequality provides the existence and finiteness of 
all its partial T derivatives of Fq(s,T, (ib) regardless to 
the sign of cr(T). Thus, using the Taylor expansion in 
powers of (T — T cep ) at any point of the interval T = T cep 



and /is < Ms(2cep), one can calculate Sh{T,/j,b) for the 
values of T > T cep which are inside the convergency ra- 
dius of the Taylor expansion. 

The other situation is for /xb > [i B (T cep ) and T > T cepi 
namely in this case above the deconfinement PT there 
must exist a weaker PT induced by the disappearance 
of the reduced surface tension coefficient. To demon- 
strate this we have solve Eq. ([5]) in the limit, when 



Indeed, using (JTUJ) , (fTTj) and (fT^J), one can write 

OA _ G 2 + u/C r _^(A,7 2 )2 7 7 / _ A 2 ~ T G 2 
dT ~ i + u/C T _i(A, 7 2 ) * uK. T -i{l,v) 
2 77 ' A 1 "" v) - AC T _i- x (l, i/)] 

(r — 1 — x) /Cr_i(l, i/) 



(25) 



where the prime denotes the partial T derivative. Note 



T approaches the curve T = T cep from above, i.e. for that the function G2 = F' + u'/C T (A, 7 2 



T 



0, and study the behavior of T deriva- 



tives of the solution of Eq. ([5]) s* for fixed values of 
fiB- For this purpose we have to evaluate the integrals 
/C T (A,7 2 ) introduced in Eq. (fT2|l . Here the notations 
A = s* - sq{T^b) and 7 2 = ~a{T) > are introduced 
for convenience. 

To avoid the unpleasant behavior for r < 2 it is con- 
venient to transform (|f 2[) further on by integrating by 
parts: 



/C T (A, 7 2 ) =9r(V )- 



(r-1) 



/C r -i(A, 7 2 )- 



7 



/C r -.(A, 7 2 ), 



(20) 



where the regular function g T (Vo) is defined as 



9r(Vo) = -, — — exp [-AV + 7 2 V] . (21) 

( T - l ) v o 

For r — a > 1 one can change the variable of integration 
v — > z/A and rewrite /C T _ a (A,7 2 ) as 



exp 



2 



/C T _ a (A, 7 2 ) = £J-a-i, dz 

V A 

& r - a - 1 K T -a(W&-") ■ (22) 

This result shows that in the limit 7^0, when the right- 
most singularity must approach sq(T,/j,b) from above, 
i.e. A — > + , the function ([2"2")) behaves as /C r _ Q (A, 7 2 ) ~ 
A T ^ a ^ 1 + 0(A T_a ). This is so because for 7 -> the ra- 
tio y 1 A~* cannot go to infinity, otherwise the function 
/C T _i (l, 7 2 A~^) , which enters into the right hand side of 
(|20[) . would diverge exponentially and this makes impos- 
sible an existence of the solution of Eq. |(SJ) for T = T cep . 
The analysis shows that for 7 — > there exist two possi- 
bilities: either 1/ = 7 2 A _>< — > Const or v = 7 2 A~^ — > 0. 
The most straightforward way to analyze these possibil- 
ities for 7 — > is to assume the following behavior 



A = A 1 a + 0{ 1 a+1 ) 
dT 8T 11 



0{l a )] 



{2k+l)Aaj a 
2 (T T ce p) 



(23) 
(24) 



and find out the a value by equating ([24]) with the T 
derivative (HI 



jq can vanish 

for a few values of \±b only. In the last step of deriving 
([23)) we used the identities ([2TJ]) and (f2"2"|) and dropped the 
non-singular terms. As we discussed above, in the limit 
7 — > the function ^ either remains a constant or van- 
ishes, then the term i/x/C t _2j«(1, v) in ([23|) is either of 
the same order as the constant /C T _x-> e (l, v) or vanishes. 
Thus, to reveal the behavior of ([25)) for 7 — > it is suf- 
ficient to find a leading term out of A 2 ~ r and 77'A 1 ^ 
and compare it with the assumption (|23[) . 

The analysis shows that for A 2_r < 77'A 1- ^ the last 
term in the right hand side of ([25]) is the leading one. 
Consequently, equating the powers of 7 of the leading 
terms in (|2"41 and ([2"5"|) , one finds 



7 



a ~ 2 ~ A 1 "* ax=2 for r < 1 



2fc + 1 



(26) 



where the last inequality follows from the fact that the 
term 77' A 1 ~* in ([25]) is the dominant one. 

Similarly, for A 2 ~ r > 77'A 1 ^ one obtains 7 Q_1 7' ~ 
A 2 ~ r and, consequently, 



for t > 1 



(r - l)(2fc + l) "" ' ~ ~ ' 2fc + 1 
Summarizing our results for 7 — > as 



(27) 



5A 



T cep 7° _ 



T-T~ 



T—T 



,T < 1 + 



2fe+l 



(28) 



r > 1 



2fc+l ' 



we can also write the expression for the second derivative 
of A as 

2fc + l _„ 

T-T ce 1 



d 2 A 
5T 2 



T-T„ 



,T< 1 



T > 1 



2k+l 



2k+l 



(29) 



The last result shows us that, depending on h and k 
values, the second derivatives of s* and sq(T,ij,b) can 
differ from each other for | < r < 2 or can be equal for 
1 < t < |. In other words, we found that at the line 



T 



T cep there exists the 2 nd order PT for § < r < 2 
and the higher order PT for 1 < r < I, which separates 
the pure QGP phase from the region of a cross-over, i.e. 
the mixed states of hadronic and QGP bags. Since it 
exists at the line of a zero surface tension, this PT will 
be called the surface induced PT. 
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FIG. 5: A schematic picture of the deconfinement phase tran- 
sition diagram (full curve) in the plane of baryonic chemical 
potential [ib and T for the 2 nd order PT at the tricritical 
endpoint (CEP). The model predicts an existence of the sur- 
face induced PT of the 2 nd or higher order (depending on 
the model parameters). This PT starts at the CEP and goes 
to higher values of T and/or fin- Here it is shown by the 
dashed curve CEP- A, if the phase diagram is endless, or by 
the dashed-dot curve CEP-B, if the phase diagram ends at 
T = 0. Below (above) each of the dashed curves the reduced 
surface tension coefficient is positive (negative), 



Since the analysis performed in the present section did 
not include any [Lb derivatives of A, it remains valid for 
the \xb dependence of the reduced surface tension coef- 
ficient, i.e. for T cep (fiB)- Only it is necessary to make a 
few comments on a possible location of the surface ten- 
sion null line T cep (/is). In principle, such a null line 
can be located anywhere, if its location does not contra- 
dict to the sufficient conditions (jTTJ) and (fl"8"|) of the 1 st 
deconfinement PT existence. Thus, the surface tension 
null line must cross the deconfinement line in the [Ib — T 
plane at a single point which is the tricritical endpoint 
(/i^ p ; T cep (/i^ p )), whereas for \ib > Ms P the null line 
should have higher temperature for the same hb than 
the deconfinement one, i.e. T , C ep(/ i s) > T c (/j,b) (see Fig 
5). Clearly, there exist two distinct cases for the surface 
tension null line: either it is endless, or it ends at zero 
temperature. But recalling that at low temperatures and 
high values of the baryonic chemical potential there may 
exist the Color-Flavor-Locked (CFL) phase [4(| , it is pos- 
sible that the null line may also cross the boundary of the 
CFL phase and, perhaps, it may create another special 
point at this intersection. 

One may wonder why this surface induced PT was not 
observed so far. The main reason is that the lattice QCD 
calculations at non-zero [1b are very difficult, and be- 
cause of this the identification of the precise location of 
the critical endpoint is highly nontrivial task [H, [i| [l(| ■ 
Therefore, the identification of the 2 nd or higher order 
PT which might be located in the vicinity of the decon- 



finement PT could be a real challenge. In addition, the 
surface induced PT may lie so close to the deconfinement 
line that it would be extremely difficult to observe it at 
the present lattices. 

VI. CONCLUSIONS AND PERSPECTIVES 

Here we suggested an analytically solvable statistical 
model which simultaneously describes the 1 st and 2 nd 
order PTs with a cross-over. The approach is general 
and can be used for more complicated parameterizations 
of the hadronic mass-volume spectrum, if in the vicinity 
of the deconfinement PT region the discrete and con- 
tinuous parts of this spectrum can be expressed in the 
form of Eqs. (TlT| and ([15)) , respectively. Also the actual 
parameterization of the QGP pressure p — Tsq(T, hb) 
was not used so far, which means that our result can be 
extended to more complicated functions, that can con- 
tain other phase transformations (chiral PT, or the PT 
to color superconducting phase) provided that the suffi- 
cient conditions (|17[) and (TT5)) for the deconfinement PT 
existence are satisfied. 

In this model the desired properties of the decon- 
finement phase diagram are achieved by accounting for 
the temperature dependent surface tension of the quark- 
gluon bags. As we showed, it is crucial for the cross-over 
existence that at T — T cep the reduced surface tension co- 
efficient vanishes and remains negative for temperatures 
above T cep . Then the deconfinement fig — T phase dia- 
gram has the 1 st PT at \i B > n° B (T cev ) for | < t < 2 , 
which degenerates into the 2 nd order PT (or higher order 
PT for | > t > 1) at hb = n B {T cev ), and a cross-over 
for < Ms < l l B(Tcep)- These two ingredients drastically 
change the critical properties of the GBM Q and resolve 
the long standing problem of a unified description of the 
I s * and 2 nd order PTs and a cross-over, which, despite 
all claims, was not resolved in Ref. Q. In addition, we 
found that at the null line of the surface tension there 
must exist the surface induced PT of the 2 ,id or higher 
order, which separates the pure QGP from the mixed 
states of hadrons and QGP bags, that coexist above the 
cross-over region (see Fig. 5). Thus, the QGBST model 
predicts that the QCD critical endpoint is the tricritical 
endpoint. It would be interesting to verify this prediction 
with the help of the lattice QCD analysis. For this one 
will need to study the behavior of the bulk and surface 
contributions to the free energy of the QGP bags and/or 
the string connecting the static quark-antiquark pair. 

In contrast to popular mean-field models the PT mech- 
anism in the present model is clear: it happens due to the 
competition of the rightmost singularities of the isobaric 
partition function. Since the GCP function of the QG- 
BST model does not depend on any (baryonic or entropy 
or energy) density, but depends exclusively on T, (1b and 
V, its phase diagram does not contain any back bend- 
ing and/or spinodal instabilities [4l| which are typical 
for the mean-field (= classical) models. The found exact 
analytical solution does not require a complicated and 



9 



artificial procedure of conjugating the two parts of the 
equation of state in the vicinity of the critical endpoint 
like it is done by hands in Refs. [H,|43[ because all this 
is automatically included in the statistical description. 

Also in the QGBST model the pressure of the de- 
confined phase is generated by the infinite bag, whereas 
the discrete part of the mass- volume spectrum plays an 
auxiliary role even above the cross-over region. There- 
fore, there is no reason to believe that any quantitative 
changes of the properties of low lying hadronic states gen- 
erated by the surrounding media (like the mass shift of 
the lo and p mesons [44j |) would be the robust signals of 
the deconfmement PT. On the other hand, the QGP bags 
created in the experiments have finite mass and volume 
and, hence, the strong discontinuities which are typical 
for the I s * order PT should be smeared out which would 
make them hardly distinguishable from the cross-over. 
Thus, to seriously discuss the signals of the I s * order 
deconfmement PT and/or the tricritical endpoint, one 
needs to solve the finite volume version of the QGBST 
model like it was done for the SMM [23] and the GBM 
[25l | . This, however, is not sufficient because, in order to 
make any reliable prediction for experiments, the finite 
volume equation of state must be used in hydrodynamic 
equations which, unfortunately, are not suited for such a 
purpose. Thus, we are facing a necessity to return to the 
foundations of heavy ion phenomenology and to modify 
them according to the requirements of the experiments. 
The present model can be considered as the next step in 
this direction. 

Although the present model has a great advantage 
compared to other models because, in principle, it can 



be formulated on the basis of the experimental data on 
the degeneracies, masses and eigen volumes of hadronic 
resonances in the spirit of Ref. ,32|, a lot of additional 
work is necessary to properly study the issues addressed 
in [45[ . Thus, above the surface tension null line the 
hadrons can coexist with QGP at high temperatures. 
Consequently, the nonrelativistic consideration of hard 
core repulsion in the present model should be modified 
to its relativistic treatment for light hadrons like it is 
suggested in [Sfl E3] . This can lead to some new effects 
discussed recently in [47| ■ Also, the realistic equation of 
state requires the inclusion of the temperature and mass 
dependent width of heavy resonances into a continuous 
part of the mass- volume spectrum which may essentially 
modify our understanding of the cross-over mechanism 

Finally, a precise temperature dependence of the sur- 
face tension coefficient of the bags should be investigated 
and its relation to the interquark string tension should 
be studied in detail. For this it will be necessary to mod- 
ify the Hills and Dales Model [13, HH m order to in- 
clude the surface deformations with the base of arbitrary 
size whereas its present formulation is suited for discrete 
clusters and, hence, for discrete bases of surface defor- 
mations. 
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